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Abstract 



Data collected at the Z resonance using the DELPHI detector at LEP are used to 
determine the charged hadron multiplicity in gluon and quark jets as a function 
of a transverse momentum-like scale. The colour factor ratio, Ca/Cp, is directly 
observed in the increase of multiplicities with that scale. The smaller than 
expected multiplicity ratio in gluon to quark jets is understood by differences 
in the hadronization of the leading quark or gluon. From the dependence of the 
charged hadron multiplicity on the opening angle in symmetric three-jet events 
the colour factor ratio is measured to be: 

^ = 2.246 ± 0.062 {stat.) ± 0.080 {syst.) ± 0.095 {theo.) 
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1 Introduction 

The gauge symmetry underlying the Lagrangian of an interaction directly determines 
the relative coupling of the vertices of the participating elementary fields. A comparison of 
the properties of quark and gluon jets, which are linked to the quark and gluon couplings, 
therefore implies a direct and intuitive test of Quantum Chromodynamics, QCD, the 
gauge theory of the strong interaction. 

Hadron production can be described via a so-called parton shower, a chain of succes- 
sive bremsstrahlung processes, followed by hadron formation which cannot be described 
perturbatively. As bremsstrahlung is directly proportional to the coupling of the radiated 
vector boson to the radiator, the ratio of the radiated gluon multiplicity from a gluon 
and quark source is expected to be asymptotically equal to the ratio of the QCD colour 
factors: Ca/Cf = 9/4 As the radiated gluons give rise to the production of hadrons, 
the increased radiation from gluons should be reflected in a higher hadron multiplicity 
and also in a stronger scaling violation of the gluon fragmentation function 

It was however noted already in the first paper comparing the multiplicities from 
gluons and quarks that this prediction does not immediately apply to the observed 
charged hadron multiplicities at finite energy as this is also influenced by differences of 
the fragmentation of the primary quark or gluon. These differences must be present 
because quarks are valence particles of the hadrons whereas gluons are not. This is 
most clearly evident from the behaviour of the gluon fragmentation function to charged 
hadrons at large scaled momentum where it is suppressed by about one order of magnitude 
compared to the quark fragmentation function . This suppression also causes a higher 
multiphcity to be expected from very low energy quark jets compared to gluon jets. 
Moreover, as low momentum, large wavelength gluons cannot resolve a hard radiated 
gluon from the initial quark-antiquark pair in the early phase of an event, soft radiation 
and correspondingly the production of low energy hadrons is further suppressed 0,^,^ 
compared to the naive expectation. In a previous publication 0] it has been shown that 
a reduction of the primary splittings of gluons compared to the perturbative expectation 
is indeed responsible for the observed small hadron multiplicity ratio between gluon and 
quark jets. 

If heavy quark jets are also included in the comparison, a further reduction of the 
multiplicity ratio is evident due to the high number of particles from the decays of the 
primary heavy particles. 

Furthermore, the definition of quark and gluon jets in three-jet events in e^e~ anni- 
hilation uses jet algorithms which combine hadrons to make jets. Low energy particles 
at large angles with respect to the original parton direction are likely to be assigned to a 
different jet. As gluon jets are initially wider than quark jets this presumably leads to a 
loss of multiplicity for gluon jets and a corresponding gain for quark jets. 

The effects discussed lead to a ratio between the charged hadron multiplicities from 
gluon and quark jets being smaller than the ratio between gluon radiation from gluons and 
from quarks. So far these effects have mainly been neglected in experimental and more 
elaborate theoretical investigations 0| . However, as we will show in this paper, at current 
energies these non-perturbative effects are still important and need to be considered in 
a proper test of the prediction |]T| that the radiated gluon-to-quark multiplicity ratio is 
equal to the colour factor ratio. 

The stronger radiation from gluons is expected to become directly evident from a 
stronger increase of the gluon jet multiplicity with the relevant energy scale as compared 
to quark jets. In this way the size of the non-perturbative terms can also be directly 
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estimated from the quark and gluon jet multiplicity at very small scales. A scale depen- 
dence of quark and gluon properties was first demonstrated in with the jet energy as 
the intuitive scale. This result was later confirmed by other measurements and 



has recently been extended to a transverse momentum-like scale |T2 



A study of the total charged multiplicity of symmetric three-jet events as function 
of the internal scales of the event avoids some of the complications mentioned above. 
A novel precision measurement of the colour factor ratio Ca/Cf can be performed by 
combining these data with a Modified Leading Log Approximation (MLLA) prediction 
of the three-jet event multiplicity [l^ which includes coherence of soft gluon radiation. 

This letter is based on a data analysis which is similar to that presented in previous 
papers We therefore have restricted the experimental discussion in section ^ to the 
relevant differences with respect to these papers. More detailed information can also be 
found in |I^ , p!5| . In section the ratio of the slopes of the mean hadron multiplicities 
in gluon and quark jets with scale is shown to be determined by the colour factor ratio 
Ca/Cf- In order to describe the data with the perturbative QCD expectations it is 
necessary to introduce additional non-perturbative offsets. This analysis is intended to 
be mainly qualitative and in many aspects it is similar to previous analyses. Then in 
section ^]2| a precision measurement of the colour factor ratio from symmetric three-jet 
events is discussed and an estimate for the difference of non-perturbative contributions 
to the quark and gluon jet multiplicity is given. Finally we summarize and conclude. 



2 Data Analysis 

The analysis presented in this letter uses the full hadronic data set collected with the 



DELPHI detector (described in [16|) at Z energies in the years 1992 to 1995. The cuts 
applied to charged and neutral particles and to events in order to select hadronic Z decays 
are identical to those given in for the qqg analysis and to |^ for the qq^ analysis. For 
the comparison of gluon and quark jets, three-jet events are clustered using the Durham 



algorithm |]T^. In addition it was required that the angles, ^2,3, between the low-energy 
jets and the leading jet are in the range from 100° to 170° (see Fig. |l]a)). Within this 
sample, events are called symmetric if 62 and 9^ are equal within some analysis-dependent 
tolerance. The leading jet is not used in the gluon or quark jet analysis. 

The identification of gluon jets by anti-tagging of heavy quark jets is identical to that 
described in Quark jets are taken from qqg events which have been depleted in 

b-quark events using an impact parameter technique. In order to achieve multiplicities 
of pure quark and gluon jet samples, the data have been corrected using purities from 
simulated events generated with JETSET 7.3 |18|] with parameters set as given in |19[. 



This is justified by the good agreement between data and simulation. Furthermore the 
model independent techniques described in 0] for symmetric events (see Fig. |l|b)) give 
results largely compatible with those obtained with the simulation correction |T^. The 
effects of the finite resolution and acceptance of the detector and of the cuts applied are 
corrected for by using a full simulation of the DELPHI detector |]TB . 



The correction for the remaining b-quark events in the qqg sample does not influence 
the slope of the measured multiplicity with scale, but only leads to a shift of its absolute 
value. 

In the simulation, quark and gluon jets are identified at "parton level" . The partons 
entering the fragmentation of a three-jet event are clustered into three jets using the 
Durham algorithm. Then for each parton jet, the number of quarks and antiquarks are 
summed where primary quarks contribute with weight +1 and antiquarks with the weight 
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Jet 2 



Asymmetric events Symmetric (Y) events 

Figure 1: Definition of event topologies and angles used throughout this analysis. The 
length of the jet lines indicates the energies. In the symmetric (Y) events (see Fig. |l]b)) 

— 1. Other quarks and gluons are assigned the weight 0. These sums are expected to 
yield +1 for quarks jets, —1 for anti-quark jets and for gluon jets. The small amount 
of events not showing this expected pattern of (+1,-1,0) was discarded. Finally, the 
parton jets were mapped to the jets at the hadron level by requiring the sum of angles 
between the parton and hadron jets to be minimal. Events exceeding a maximum angle 
between the parton and jet directions were also rejected. At large opening angles the 
influence of these rejections is found to be about 3% increasing at low opening angles. 

The gluon jet purities vary from 95% for low energy gluons to 46% for the highest 
energy gluons. The few bins with lower purities have been excluded from the analysis. 
The quark purities range from 43% to 81%. 

For the analysis of the multiplicity of symmetric three-jet events, all events were forced 
to three jets using the Durham algorithm without a minimal ycut- The angles between the 
jets were then used to rescale the jet momenta to the centre-of-mass energy as described in 
Ip. Symmetric events were selected by demanding that 62 be equal to 6^ within 2e. Here 
e is half the angular bin width of 9i taken to be 3°. The analysis has been performed 
for events of all flavours as well as for b-depleted events. In both cases the measured 
multiplicity was corrected for track losses due to detector effects and cuts applied. The 
correction factor was calculated as ratio of generated over accepted multiplicity using 
simulated events. It varies smoothly, from 1.25 at small 9i to 1.32 at large 9i. 

3 Results 

3.1 Comparison of Multiplicities in Gluon and Quark Jets 

In order to determine a scale dependence, the scale underlying the physics process 
needs to be specified. The actual physical scale is necessarily proportional to any variation 
of an outer scale like the centre-of-mass energy. As usually only the relative change in 
scale matters, this outer scale can therefore be used instead of the physical scale. For 
this analysis the situation is different. The jets entering the analysis stem from Z decays 
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and thus from a fixed centre-of-mass energy. So the relevant scales have to be determined 
from the properties of the jets and the event topology. From the above discussion the 
scale has to be proportional to the jet energy because this quantity scales with the energy 
in the centre-of-mass system for similar events. Studies of hadron production in processes 
with non-trivial topology have shown that the characteristics of the parton cascade prove 
to depend mainly on the hardness of the process producing the jet p.p(]|: 

e 

K = Ejetsm- . (1) 

Ejet is the energy of the jet and 6 its angle to the closest jet. This scale definition corre- 
sponds to the beam energy in two-jet events. It is similar to the transverse momentum 
of the jet and also related to y/i/cut as used by the jet algorithms. It is also used as the 
scale in the calculation of the energy dependence of the hadron multiplicity in e'^e~ an- 



nihilation [0j2^ to take into account the leading effect of coherence. It should, however, 
be noted that several scales may be relevant in multi-jet events. Hence using k is an 
approximation. A similar scale, namely the geometric mean of the scales of the gluon jet 
with respect to both quark jets while using Eqn. |l] for the quark jets, has recently been 
used in a study of quark and gluon jet multiplicities [|12 . 



As stated in the introduction we want to gain information on the relative colour charges 
of quarks and gluons from the rate of change of the multiplicities with scale. Assuming 
the validity of the perturbative QCD prediction, the ratio of the charged multiplicities 
of gluon and quark jets, Ngiuon/Nquarki has to approach a constant value (approximately 
the colour factor ratio) at large scale. This trivially implies that the ratio of the slopes 
of quark and gluon jet multiplicities also approaches the same limit. This fact is a direct 
consequence of de rHopital's rule ||2^ and is also directly evident from the linearity of 
the derivative: 

at large scale: Ngiuon{i^) = C ■ Nquark{i^) dNgiuon/ d^. ^ (^2) 

0'-'^ quark / (^1^ 

i.e. the QCD prediction for the ratio of multiplicities applies equally well to the ratio 
of the slopes of the multiplicities. In fact it is to be expected that the slope ratio is 
closer to the QCD prediction than the multiphcity ratio as it should be less affected by 
non-perturbative effects. 

This effect has been cross-checked using the HERWIG model which allows the 



number of colours to be changed and thus by SU(n) group relations, the colour factor 
ratio Ca/Cf- The predictions of HERWIG are found to follow directly the expectation 
of the right hand side of Eqn. This has also been confirmed in a recent theoretical 
calculation of this quantity in the framework of the dipole model. 

Fig. ^) shows the multiplicity in quark and gluon jets as a function of the hardness 
scale K. For both multiplicities an approximately logarithmic increase with k is observed 
which is about twice as big for gluon jets as for quark jets, thus already strikingly con- 
firming the QCD prediction. 

A stronger increase of the gluon jet multiplicity was already noted in a previous paper 
1^, where the jet energy was chosen as scale. Meanwhile this observation has been 
confirmed also by other measurements P, pil|JTT| and has been extended to different scales 
|T^. Fragmentation models (not shown) predict an increase of the multiplicities which is 



in good agreement with the data. 

In order to obtain quantitative information from the data shown in Fig. |^a), the 
following ansatz was fitted to the data: 
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Figure 2: a) Average charged particle multiplicity for light quark and gluon jets as 
function of k fitted with Eqn. ^; b) ratio of the gluon to quark jet multiplicity; the full 
line shows the ratio of the functions fitted to the data in a), the dashed curve is the 
ratio of the slopes of the fits in a). All curves are extrapolated to the edges of the plot 
by the dotted lines. Also included are measurements of the multiplicity ratio of some 
other experiments |iyjTl|. The grey band shown with the slope ratio indicates the error 



estimated by varying all fit parameters within their errors. 



<Ng>{K) = + Np,rt{^) ■ r{K) (3) 

Here Nq^ are non-perturbative terms introduced to account for the differences in the 
fragmentation of the leading quark or gluon as discussed in detail in the introduction. 
These terms are assumed to be constant. Npert is the perturbative prediction for the 



hadron multiplicity as given in ||2T[| : 



Np,rt{K) = K ■ {as{K)f ■ exp | ] ■ [1 + 0{^s)] (4) 

'^-A^3joV-C-J ^ '^^^ ^ /^o-ll--/ . 

A first and a second order as have been used with this expression with the number of 
active flavours, nj, equal to five. An alternative prediction has been given in using 
the limited spectrum approach: 

iVpert = K • r(s) (1)' ""/i+bW (5) 

^ 33 + 2/9n/ , 
^= 33-2n, ' ^ = l°gA^7o ; 7o 
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Here a first order has always been used with Uf taken as three T is the Gamma- 



function and Ib the modified Bessel-function. K is a non-perturbative scale factor. The 
QCD scale parameter A enters into the definition of q;s(k^/A^) |Q. The numerical values 
of K and A are not expected to be the same in Eqns. ^ and | as different approximations 
are used. Finally: 

C 

r{K) = 7^(1 - ri7o - r27o) (6) 



with: 




2nfCF\ ri /25 3 Uf njC 



- I ^ + - j ' 6 1^ 8 ACa Cl ) 

is the perturbative prediction for the multiplicity ratio in back-to-back gluon to back- 
to-back quark jets. The terms proportional to ri (r2) correspond to the NLO (NNLO) 
prediction. Numerically they correspond to corrections of about 8% and 1% respectively. 
The smallness of the higher order corrections indicates that the perturbative series of the 
gluon-to-quark multiplicity ratio converges rapidly. 

The fits represent the data well. The fit range has not been extended to too small 
scales as here a contribution of initial two-jet events might bias the multiplicities to lower 
values. Parameters of the fits for this specific choice of scale and jet selection are given in 
Tab. |T]. No estimate of systematic error is given as this analysis is intended to be mainly 
qualitative. The fit parameters should not be compared directly to those parameters 
usually obtained from overall events in e^e~ annihilation. The normalization factor, 
differs strongly due to the differences in the multiplicity in jets and overall events. 
Furthermore, the introduction of non-perturbative offsets leads to a strong reduction of 
the values of the effective scale parameter A. This is also observed if the e^e~ multiplicity 
is fitted including an offset term, which could be reasonable in this case also. 



Parameter 


iVpert from 
Eqn. 1 


Npert from 
Eqn. 1 


A[GeV] 


0.032 ± 0.011 


0.011 ± 0.004 


K 


0.005 ± 0.001 


0.12 ± 0.02 


Ca/ Cf 


2.12 ± 0.10 


2.15 ± 0.10 


K 


2.82 ± 0.14 


3.12 ± 0.20 




0.73 ± 0.21 


1.43 ± 0.31 


xVn-d.f. 


0.61 


0.65 



Table 1: Results of the fits of the quark and gluon jet multiplicities as a function of k. 



Using an identical scale definition for quark and gluon jets also allows the gluon-to- 
quark jet multiplicity ratio to be directly evaluated as function of this scale. Fig. ||b) 
shows this ratio as calculated from data and the fits as function of the hardness scale as 
well as the ratio of the slopes of the fits. The ratio of the multiplicities increases from 
about 1.15 at small scale to about 1.4 at the highest scales measured. The measurement 
1^ performed in T(IS') — > -^gg decays at small sea iJ, and of "inclusive" gluons |lT 



at large scale, agree quite well with the expectation from the fits. The corresponding 



hardness scale for the data at the highest scale [O] has been estimated from the average 



^Half of the gg invariant mass is taken as the equivalent scale. 
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gluon energy and the angle cuts given in fTTI]. The good agreement of the "inclusive' 



gluon measurement also implies that angular ordering effects are relevant in this case. 

The ratio of the slopes for the different fits is almost 2 corresponding to a colour factor 
ratio of Ca/Cf = 2.12 ± 0.10, well compatible with the QCD expectation. 

The fits further indicate that for very small scale the multiplicity of quark jets is bigger 
than that of gluon jets. Consequently the constant terms contributing to the multiplicity 
due to the primary gluon or quark fragmentation are larger for quarks (see Tab. |I|). The 
difference of these terms is about 2. Taking the scale choice made in |]T^ leads to about 
a 20% increase of the measured colour factor ratio and a corresponding increase in the 
difference of the non-perturbative constants to 4.2. 

It is instructive here to estimate a lower limit for the difference of the non-perturbative 
terms from the behaviour of the gluon and quark fragmentation functions . Due to lead- 
ing particle effects the fragmentation function of the quark outreaches the fragmentation 
function of the gluon at high values of xe- Taking the shape of the gluon fragmentation 
function as unbiased by the leading particle effect and assuming the overall multiplicity of 
gluon jets roughly as twice as big as of quark jets, one gets an estimate for the lower limit 
of additional multiplicity in quark jets by integrating the difference between the quark 
and the halved gluon fragmentation function in the X£;-region where the fragmentation 
function of the gluon is below that of the quark. This yields Nq — Nq > 0.61 ± 0.02 from 
Y and Nq — A^'q > 0.58 ± 0.05 from so-called Mercedes events 0. It should be noted 
here, that the leading particle effect still influences the multiplicity at even lower scaled 
hadron energies. The region of small hadron energy contributes most to the multiplicity. 
Therefore the estimated limit presumably is much smaller than the actual value of Nq. 

At first sight a difference of the constant terms of the order of ~2 units in charged 
multiplicity looks unexpectedly large. However, these constants also include the effects 
of the jet clustering. Furthermore, stable hadron production to a large extent proceeds 
via resonance decays, so that the observed difference may only correspond to a difference 
of about one primary particle. The larger constant term for quarks compared to gluons 
explains the different behaviour of the ratio of multiplicities and the slope ratio in Fig. ^d). 

The observed behaviour would be expected from non-perturbative effects of the frag- 
mentation in the leading quark or gluon. In the cluster fragmentation model, an ad- 
ditional gluon to quark-antiquark splitting is needed in the fragmentation of a gluon 
compared to that of a quark. 



3.2 Precise Determination of Ca/Cf from Multiplicities in 
Three-Jet Events 

The analysis presented so far, as in most other comparisons of quark and gluon jet 
multiplicities, has the disadvantage of relying on the association of (maybe low energy) 
particles to jets. Clearly this involves severe ambiguities and specifically does not consider 
coherent soft gluon radiation from the initial qqg ensemble. This can be avoided and a 
precise measurement can be obtained by studying the dependence of the total charged 
multiplicity in three-jet events as function of the quark and gluon scales. In fact there is 
a definite MLLA prediction [|13| for this multiplicity Nggg-. 

+ (7) 



A^, 



ggg 



2Ng{Y;g) + NgiY* 
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with the scale variables: 



Y* — In 



5 



In 



{PgPg)iPgPg) _ ^ 
^ 2A2(p,p,) 2A 



(8) 



Nq(Y*g) and A'gCi^*) describe the scale dependence of the multiplicity for quark or gluon 
jets, respectively. A is a scale parameter and the Pq,q,g are the four-momenta of the 
quarks and the gluon. The three-jet multiplicity depends on the quark energy, E*, in the 
centre-of-mass system of the quark-antiquark pair and on the transverse momentum scale 
of the gluon, pj;. For comparison with data, this is expressed in as a dependence 
on the measured multiplicity in e+e~ events, N^+e-, and the colour factor ratio as given 
in Eqn. |^. In addition, we again choose to add a constant term, A^^o, to account for 
differences in the fragmentation of quarks and gluons as discussed above. Thus, omitting 
correction terms: 



N, 



ggg 



(9) 



Although at first sight this appears to be the incoherent sum of the multiplicity of the two 
quark jets and the gluon jet, this formula includes coherence effects in the exact definition 
of the scales of the N^+e- terms [P7|. Nevertheless, subtracting the non-perturbative 
term No within the curly brackets gives a physical interpretation for Nq as the additional 
multiplicity in quark jets due to the leading particle effect, which is contained in the 
measured Nf>+e- and has to be subtracted to get the gluon contribution to the multiplicity. 

In principle Eqn. ^ still requires the determination of the quark-antiquark and gluon 
scales independently. However, in symmetric Y-type events (see Fig. [l|b)) both scales can 
be expressed as functions of the opening angle 6i only by initially assuming that the gluon 
jet is not the most energetic one. E*^ oc EgEgSiii^ 9^/2 for this type of event is almost 
constant (see upper full curve in Fig. |^)) at fixed centre-of-mass energy. However, pj-, 
increases approximately linearly with the opening angle as it is proportional to the gluon 
transverse momentum. As the multiplicity change corresponding to the change of E* 
corresponds only to about —2, the 9 dependence of the three-jet multiplicity therefore 
mainly measures the scale dependence of the multiplicity of the gluon jet. 

In a fraction of the events (which strongly increases with opening angle) the gluon jet 
is the most energetic jet. This can be corrected for in different ways when fitting Eqn. ^ to 
the data using Monte Carlo simulation. Assuming an approximately logarithmic increase 
of the multiplicity with scale, which is well supported by the data, the average scale at 
a given opening angle can be expressed as the geometric mean of the cases where the 
gluon initiates the most energetic jet and where it does not. These corrected scales are 
shown as the points in Fig. ^). The correction first increases with the opening angle but 
then decreases again and vanishes for fully symmetric events. Alternatively, the fraction 
of events when the gluon initiates the most energetic jet can be considered separately in 
Eqn. |. 

To obtain information on the colour factor ratio Ca/Cf, the scale dependence of the 
three-jet multiplicity has to be compared to the multiplicity in all e~^e~ events. This has 
been chosen to be taken from the DELPHI measurements with hard photon radiation for 
energies below the Z mass and at 184 GeV PDI and the LEP combined measurements 
at the intermediate energies [^. For studies of systematic errors, data from lower en- 
ergy e~^e~ experiments have also been used. The DELPHI multiplicities in events 
with hard photon radiation have been extracted as described in |]8|JT^, but using the 
full statistics now available. Small energy dependent corrections (2 — 4%) to the e~^e~ 
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Figure 3: a) Variation of the scales 2E* and as function of the opening angle 9i in 
symmetric three-jet events. The functions are the analytic expectation. The points 
include a correction (calculated with JETSET 7.3) for the cases where the gluon forms 
the most energetic jet. The lines matching the points are polynomials fitted to obtain 
continuous values. 

b) Charged hadron multiplicity as a function of the centre-of-mass energy of the gg-pair 
fitted with the perturbative predictions Eqs. ^ or ^. 

c) Charged hadron multiplicity in symmetric three-jet events as a function of the opening 
angle. The dashed curve is the prediction using the ansatz Eqn. || setting Ca/Cf to its 
default value and omitting the constant offset, Nq. The full curve is a fit of the full ansatz 
Eqn. H to the data treating Ca/Cf and A'o as free parameters. 

d) Stability of the result for Ca/Cf against variation of the smallest opening angle used in 
the fit as well as x^/^df and the probability of these fits. The dash-dotted horizontal 
line in the upper half shows the QCD expectation for Ca/Cf with the dotted lines 
representing variations of ±10%. 

The DELPHI data of Fig. ^ b) and c) will be made available in the Durham/RAL 

database 
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multiplicities were applied to correct for the varying contribution of b quarks. The mul- 
tiplicities obtained were fitted with the perturbative predictions, Eqs. ^or^], see Fig^). 
Both calculations describe the data equally well. The parameters of the fits are given in 
the upper part of Tab. 0. 



Parameter 


Np,rt from 121P 
(Eqn.i 


Np,rt from |2| 
(Eqn.D 


relevant 
data 


A 
K 

xVn.d.f. 


0.275 ± 0.070 
0.026 ± 0.003 
1.180 


0.061 ± 0.015 
0.606 ± 0.062 
1.183 


data from 
e~^e~ and qq'j 


Ca/ Cf 
No 

xVn.d.f. 


2.251 ± 0.063 
1.40 ± 0.10 
0.998 


2.242 ± 0.062 
1.40 ± 0.10 
1.004 


data from 
symmetric 
3 jet events 



Table 2: Result of the fits of the e^e multiplicity (upper part) and the three-jet event 
multiplicity (lower part). 

The measured, fully corrected multiplicity in all symmetric three-jet events as function 
of the opening angle is shown in Fig. ^). A strong increase of the multiplicity from values 
of around 18 for small opening angle to about 29 at opening angles of 120° (corresponding 
to fully symmetric events) is observed. Omitting the non-perturbative term, A'o, in Eqn. |^ 
and setting Ca/Cf to its expected value predicts a similar increase in multiplicity over 
this angular range (dashed curve in ^)). The prediction is however higher by about three 
units of charged multiplicity. This discrepancy is expected from the previously obtained 
result due to differences in the fragmentation of the leading quark or gluon. 

At small angles the difference between the primary QCD expectation and the mea- 
surement increases. Studies using Monte Carlo models have shown that this is mainly 
due to genuine two jet events which have been clustered as symmetric three-jet events. 
The models indicate that this contribution becomes small for angles above 30°. 

Fitting the full ansatz |] to the three-jet multiplicity data at angles 9 > 30°, using 
the two parameterizations in Eqs. |^ and ^ of the multiplicity in e"'"e~ events with their 
parameters fixed as given in Tab. ^ but varying Ca/Cf and A^o, yields: 

2.251 ±0.063 (10) 

2.242 ±0.062. (11) 

The result confirms with great precision the QCD expectation [|l| that the ratio of the 
radiated multiplicity from gluon and quark jets is given by the colour factor ratio Ca/Cf- 
This result also implies that the proportionality of the number of gluons to hadrons |jl| 
e.g. Local Hadron Parton Duality (LPHD) |^ applies extremely precisely if only the 
radiated gluons from a quark or gluon are considered. 

The offset term A'o is bigger if only b-depleted events are used. The central result for 
Ca/Cf, however, remains unchanged within errors. This is due to the fact that Ca/Cf 
is measured from the change of multiplicity in three-jet events with opening angle and 
not from the absolute multiplicity. 

The correctness of the ansatz Eqn. ^ and the bias introduced by two-jet events at 
small 6i, were further checked by varying the lowest angle used in the fit. The resulting 



Ca 
Cf 
Ca 
Cf 
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value for CaICf, the /Ndj and the probabihty of the fit are shown in Fig. ^jd). It 
is observed that for 9i > 30° satisfactory fits are obtained. For this angular range the 
fitted value of Ca/Cf is stable within errors. 

Systematic uncertainties of the above result for the colour factor ratio due to uncer- 
tainties in the three-jet multiplicity data as well as in the parameterization of the e^e~ 
charged multiplicity and in the theoretical predictions are considered. To obtain system- 
atic errors interpretable like statistical errors, half the difference in the value obtained for 
CaICp when a parameter is modified from its central value (see below) is quoted as the 
systematic uncertainty. All relative systematic errors are collected in Tab. ^ 



Source 


Sys. error 


combined 


combined 


total 


Experimental uncertainties 


±5.52% 


1. Min. particle momentum 

2. Min. angle of jet w.r.t. beam 

3. Min. number of tracks per jet 

4. Corr. for gluon in jet 1 


± 0.42 % 
± 0.38 % 
± 0.02 % 
± 0.11 % 


± 0.58% 


± 3.55 % 


5. jet algorithms 




± 1.39 % 


6. e+e" data sets 

7. Fit function 

8. binning and range of fit 


± 0.90 % 
± 0.02 % 
± 3.08 % 


± 3.21 % 


Theoretical uncertainties 


9. Variation oi rif 

10. Calculation in lst/2nd order 

11. Setting Ca fixed 


± 1.51 % 
± 3.95 % 
± 0.08 % 




± 4.23% 



Table 3: Systematic uncertainties on Ca/Cf as derived from three-jet event multiplicities 

Results for Ca/ Cf obtained from the individual data sets corresponding to the different 
years of data-taking as well as from b-depleted events were found to be fully compatible 
within the statistical error. To estimate uncertainties in the three-jet multiplicity the 
following cuts which are sensitive to misrepresentation of the data by the Monte Carlo 
simulation have been varied. 

1. Cut on the minimal particle momentum: 

the cut on the minimal particle momentum has been lowered from 400 MeV to 
200 MeV and raised to 600 MeV. 

2. Minimum angle of each jet with respect to the beam axis: 

this cut has been increased from 30° to 40° to test for a possible bias due to the 
limited angular acceptance. 

3. Minimum number of particles per jet: 

the minimum number of particles per jet has been increased from 2 to 4 in order to 
reject events which may not have a clear three-jet structure. 

4. Correction for gluon in leading jet: 

both methods of correction were compared to account for gluons in the most energetic 
jet. Furthermore the requirements for the mapping of the parton to the hadron level 
for defining the gluon jet have been varied. 

To check the stability of the result for different choices of jet algorithms the results 
obtained for a large sample of events generated with JETSET have been compared with: 
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5. Alternative jet algorithms: 

the angular ordered Durham algorithm, LUCLUS without particle reassignment, 
JADE and Geneva |^ were applied alternatively to Durham on a large statistics 



Monte Carlo sample. The results for Durham, angular ordered Durham and LU- 
CLUS agree reasonably. The spread among the results was taken as error. The 
JADE and Geneva algorithm which are known to tend to form so-called junk jets 
33 show stronger deviations. 



The following systematic uncertainties arise from uncertainties in the experimental 
input other than from the three-jet multiplicities and from choices made for the fits of 
iVg+e-. These uncertainties are considered as experimental systematic uncertainties. 

6. Input of parameterization of N^+g-^y/s): 

to estimate the influence of an uncertainty in Ng+e- , different choices of input data 
were compared: 

• DELPHI multiplicities for 184 GeV and from Z decays with hard photons com- 
bined with LEP data for 90 GeV < < 180 GeV ; 

• DELPHI multiplicities from Z decays with hard photons; 

• e~^e~ data taken at low centre-of-mass energies (TASSO, TPC, MARK-II, HRS, 
AMY); 

• all available e+e" data between 10 GeV and 184 GeV (TASSO, TPC, MARK-II, 
HRS, AMY, LEP combined, DELPHI). 

7. Choice of prediction used for fit: 

the fit functions ^ and | were used alternatively. For consistency here Uf = 5 and a 
second order as was used. 

8. Variation of the fitted range: 

the lower limit of the angular range used in the fit was varied between 24° and 36° 
as well as changing half the bin width, e, from 2.5° to 5°. 

Finally, systematic errors due to uncertainties in the theoretical prediction were con- 
sidered. 

9. Variation of n/: 

the number of active quarks, n/, [^] relevant for the hadronic final state is uncertain. 
Uf therefore has been varied from 3 to 5. 
10. Order of calculation (LO - NNLO): 

the prediction r(K) (Eqn. |^) has been calculated for back-to-back quarks or gluons. 
As the jets are well separated it is expected to apply for this analysis also. When 
the gluon recoils with respect to the quarks the prediction is exact. In addition 
coherence effects (angular ordering) are taken into account in the definition of the 
scales E* and p^. 

As the coupling for the triple-gluon vertex is bigger than the coupling of all other 
vertices it is clear that the correction will lower the gluon-to-quark multiplicity ratio 
as in the case of Eqn. ^. The validity of the correction is therefore assumed for 
the whole range of angles considered. Conservatively, half of the difference obtained 
with the lowest order prediction r = Ca/Cf and the NNLO prediction is considered 
as systematic uncertainty. A leading order was used for the lowest order prediction 
and a second order in the other case. Considering that in the three-jet events 
mainly the gluon scale is varied, the resulting error estimate agrees with that 
given in Eqn. |^. 
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11. Quantities influencing Ca/Cp: 

for tlie central result, Ca/Cp has been assumed variable in Eqn. ^ only. The sta- 
bility of the result was checked by also leaving Ca variable in some or all of the 
parameterizations of as and N^+g-. 

To check in how far the offset term A^o is constant, A^^o has been extracted for each 
6'i-bin individually fixing Ca to its default value. The individual results are consistent 
with the average value and no trend is observed. 

Alternatively to Eqn. ^, Eqn. ^ has been fitted to the data, where the 0{as) cor- 
rection factor has been parameterized as (1 + cas{pi))- This leads to the same fit re- 
sults for CaICp and Eqn. ||, which implies that both corrections r(pf) ■ A'q and 
2Nq{Y*^ + Ng{Y*) ■ casipi) well as the values obtained for Ca/Cp agree within 
±1%. It should, however, be stressed that the behaviour of the fragmentation function 
requires the presence of a non-perturbative offset term. 

The prediction of the multiplicity ratio given by has been tried as an alternative 



to Eqn. ^ Although this calculation takes recoil effects into account, a non-perturbative 
offset term is still required. The prediction differs by about 10% from in the NNLO 
term. As it does not reproduce the colour factor ratio contained in the fragmentation 
models which describe the data well, it has not been applied in this analysis. 



Averaging the results given in Eqns. and 11 and adding in quadrature the system- 



atic errors summarized in Tab. ^ gives the following final result: 
C 

— ^ = 2.246 ± 0.062 {stat.) ± 0.080 {syst.) ± 0.095 {theo.) (12) 
Cp 

This result confirms the QCD expectation that gluon bremsstrahlung is stronger from 
gluons than from quarks by the colour factor ratio Ca/Cp and is direct evidence for the 
triple-gluon coupling. 

This measurement yields the most precise result obtained so far for the colour factor 



ratio Ca/Cp- Even the best measurements from four-jet angular distributions suffer 
from the relatively small number of four-jet events available. Furthermore, many of these 
measurements specify no theoretical systematic error as they so far rely on leading order 
calculations. It is remarkable that this measurement of Ca/Cp is performed from truly 
hadronic quantities, the charged multiplicities. Jets, i.e. partonic quantities only enter 
indirectly via the definition of the scales E* and p^. 

In order to illustrate comprehensively the contents of the measurement of the three-jet 
multiplicity we compare in Fig. |^ the multiplicity corresponding to a gg and a qq final 
state. The qq multiplicity is taken to be the multiplicity measured in e^e~ annihilation 
corrected for the bb contribution as described above. The gg multiplicity at low scale 



values is taken from the CLEO measurement [[TQ], for which no systematic error was 
specified. At higher scale, twice the difference of the three-jet multiplicity and the qq 
term (the first term in Eqn. ^ is interpreted as the gg multiplicity. The gg data should be 
extendable to higher energies by measuring the multiplicity in pp scattering as a function 
of the transverse energy. The dashed curve through the qq points is a fit of the prediction 
according to Eqns. ^ or |^. The gg line is the perturbative expectation for back-to-back 
gluons according to the second term of Eqn. 0. A'o is taken from Eqn. O. In principle 



A'o is a property of the complete three-jet event, so it is unclear if the subtraction of the 
full amount of A'o is justified in order to obtain the gluon jet multiplicity. However, this 
only introduces a constant shift in the ^^gg event" multiplicity, the scale dependence of 
the gluon jet multiplicity remains unaltered. The plot shows again that the increase of 
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Figure 4: Comparison of the charged hadron multiphcity for an initial qq and a gg pair 
as function of the scale. The dashed curve is a fit according to Eqns. ^ or |, the full line 
is twice the second term of Eqn. ^. The grey band indicates the uncertainty due to the 
error of iVo. The DELPHI gg data will be made available in the Durham/RAL database 
31. 



the gg multiplicity with scale is about twice as big as in the qq case, illustrating the large 
gluon-to-quark colour factor ratio Ca/Cf- 

It is of interest to present also a dedicated measurement of the non-perturbative pa- 
rameter Nq. In order to obtain this value, b-depleted events have been used. A fit of the 
three-jet event multiplicity has then been performed with A^^o as the only free parameter. 
Ca/Cf has been set to its default value. The parameterization of the e+e~ multiplicity 
according to Eqn. || uses the low energy e^e~ data as input. The fit yields: 

A'o = 1.91 ± 0.03(stat.) ± 0.33(syst.) (13) 

The systematic error was estimated as for Ca/Cf- Furthermore a normalization error 
due to the multiplicity in e^e~ events has been added in quadrature. This error has been 
assumed to be given by the error of the precise average multiplicity at the Z resonance 
| 37| . The actual value of A^o ~ 2 corresponds to about one primary particle (see also 



section 3.1). This is indeed a reasonable value which had already been expected in |l| 
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4 Summary 

In summary, the dependence of the charged particle multiplicity in quark and gluon 
jets on the transverse momentum-like scale has been investigated and the charged hadron 
multiplicity in symmetric three-jet events has been measured as a function of the opening 
angle 9i. 

The ratio of the variations of gluon and quark jet multiplicities with scale agrees with 
the QCD expectation and directly reflects the higher colour charge of gluons compared 
to quarks. This can also be interpreted as direct evidence for the triple-gluon coupling, 
one of the basic ingredients of QCD. It is of special importance that this evidence is 
due to very soft radiated gluons and therefore complementary to the measurement of the 
triple-gluon coupling in four-jet events at large momentum transfer. 

The increase of the gluon to quark jet multiplicity ratio with increasing scale is under- 
stood as being due to a difference in the fragmentation of the leading quark or gluon. The 
simultaneous description of the quark and gluon jet multiplicities with scale also supports 
the Local Parton Hadron Duality hypothesis |Q although large non-perturbative terms 



for the leading quark or gluon are responsible for the observed relatively small gluon to 
quark jet multiplicity ratio. 

Using the novel method of measuring the evolution of the multiplicity in symmetric 
three-jet events with their opening angle, a precise result for the colour factor ratio is 
obtained: 

= 2.246 ± 0.062 (stat.) ± 0.080 [syst.) ± 0.095 (theo.) 

It is superior in precision to the best measurements from four-jet events Finally it is 
remarkable that this measurement is directly performed from truly hadronic quantities. 
Jets only enter indirectly via the definition of the energy scale of the quark-antiquark pair 
and the transverse momentum scale of the gluon. These scales are calculated directly 
from the jet angles. 
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